NAS A/TM - 2005-2 1 3 848 



A “Kane’s Dynamics” Model for the 
Active Rack Isolation System 

Part Three: Addition of Umbilicals to the Nonlinear Model 

J.K. Rupert 

Dynetics, Inc., Huntsville, Alabama 
R.D. Hampton 

United States Military Academy, West Point, New York 
G.S. Beech 

Marshall Space Flight Center, Marshall Space Flight Center, Alabama 


February 2005 


The NASA STI Program Office... in Profile 


Since its founding, NASA has been dedicated to 
the advancement of aeronautics and space 
science. The NASA Scientific and Technical 
Information (STI) Program Office plays a key 
part in helping NASA maintain this important 
role. 

The NASA STI Program Office is operated by 
Langley Research Center, the lead center for 
NASA’s scientific and technical information. The 
NASA STI Program Office provides access to 
the NASA STI Database, the largest collection of 
aeronautical and space science STI in the world. 
The Program Office is also NASA’s institutional 
mechanism for disseminating the results of its 
research and development activities. These results 
are published by NASA in the NASA STI Report 
Series, which includes the following report types: 

• TECHNICAL PUBLICATION. Reports of 
completed research or a major significant 
phase of research that present the results of 
NASA programs and include extensive data 
or theoretical analysis. Includes compilations 
of significant scientific and technical data 
and information deemed to be of continuing 
reference value. NASA’s counterpart of peer- 
reviewed formal professional papers but has less 
stringent limitations on manuscript length and 
extent of graphic presentations. 

• TECHNICAL MEMORANDUM. Scientific 
and technical findings that are preliminary or of 
specialized interest, e.g., quick release reports, 
working papers, and bibliographies that contain 
minimal annotation. Does not contain extensive 
analysis. 

• CONTRACTOR REPORT. Scientific and 
technical findings by NASA-sponsored 
contractors and grantees. 


• CONFERENCE PUBLICATION. Collected 
papers from scientific and technical conferences, 
symposia, seminars, or other meetings sponsored 
or cosponsored by NASA. 

• SPECIAL PUBLICATION. Scientific, technical, 
or historical information from NASA programs, 
projects, and mission, often concerned with 
subjects having substantial public interest. 

• TECHNICAL TRANSLATION. 
English-language translations of foreign 
scientific and technical material pertinent to 
NASA’s mission. 

Specialized services that complement the STI 
Program Office’s diverse offerings include creating 
custom thesauri, building customized databases, 
organizing and publishing research results. . .even 
providing videos. 

For more information about the NASA STI Program 
Office, see the following: 

• Access the NASA STI Program Home Page at 
http:/ lwww.sti.nasa.gov 

• E-mail your question via the Internet to 
help@sti.nasa.gov 

• Fax your question to the NASA Access Help 
Desk at 301-621-0134 

• Telephone the NASA Access Help Desk at 
301-621-0390 

• Write to: 

NASA Access Help Desk 

NASA Center for AeroSpace Information 

7121 Standard Drive 

Hanover, MD 21076-1320 

301-621-0390 


NAS A/TM - 2005-2 13848 



A “Kane’s Dynamics” Model for the 
Active Rack Isolation System 

Part Three: Addition of Umbilicals to the Nonlinear Model 

J.K. Rupert 

Dynetics, Inc., Huntsville, Alabama 
R.D. Hampton 

United States Military Academy, West Point, New York 
G.S. Beech 

Marshall Space Flight Center, Marshall Space Flight Center, Alabama 


National Aeronautics and 
Space Administration 

Marshall Space Flight Center • MSFC, Alabama 35812 


February 2005 


TRADEMARKS 


Trade names and trademarks are used in this report for identification only. This usage does not constitute an official 
endorsement, either expressed or implied, by the National Aeronautics and Space Administration. 


Available from: 


NASA Center for AeroSpace Information 
7121 Standard Drive 
Hanover, MD 21076-1320 
301-621-0390 


National Technical Information Service 
5285 Port Royal Road 
Springfield, VA 22161 
703-487^1650 


11 


TABLE OF CONTENTS 


1. INTRODUCTION 1 

2. COORDINATE SYSTEMS 3 

3. GENERALIZED COORDINATES AND GENERALIZED SPEEDS 5 

4. ANGULAR VELOCITIES OF REFERENCE FRAMES AND RIGID BODIES 6 

5. GENERALIZED ACTIVE FORCE CONTRIBUTIONS DUE TO THE UMBILICALS 8 

5.1 Generalized Active Force Equations 8 

5.2 Partial Velocities Sl v^ u J 9 

5.3 Partial Angular Velocities 5| 0) r F ' 11 

5.4 General Form for Umbilical Forces and Moments 11 

5.5 Umbilical Elongations and Elongation Rates 12 

5.6 Umbilical Angles of Twist and Twist Rates 14 

6. IMPLEMENTATION IN AUTOLEV 17 

7. MODEL VERIFICATION 18 

7.1 Onboard Impulsive-Disturbance Force, No Damping 19 

7.2 Onboard Sinusoidal-Disturbance Force, With Damping 19 

7.3 Onboard Impulsive-Disturbance Moment, No Damping 21 

7.4 Onboard Disturbance Moment, With Damping 22 

7.5 Off-Board Translational Disturbance 24 

7.6 Off-Board Rotational Disturbance 27 

7.7 Comparisons Using Nondiagonal Stiffness and Damping Matrices 29 

8. CONCLUSION 30 


REFERENCES 


31 


LIST OF FIGURES 


1. Detailed diagram of the umbilical assembly 1 

2. Kinematic diagram, including the z'th actuator assembly and the umbilical 3 

3. CAD-based technique for determining umbilical attachment locations 17 

4. Translational single-degree-of-freedom truth model 18 

5. Translational displacement due to onboard impulsive-disturbance force 19 

6. Translational displacement due to onboard sinusoidal-disturbance force 

with = 1 and £ = 1 20 

co n 

7. Translational displacement due to onboard sinusoidal-disturbance force 

with - 5 and £ = 1 20 

co n 

8. Translational displacement due to onboard sinusoidal-disturbance force 

with = 1 and £ = 0.25 21 

(On 

9. Rotational displacement due to onboard impulsive-disturbance moment 22 

10. Rotational displacement due to onboard sinusoidal-disturbance moment 

with = 1 and £ = 1 23 

(O n 

11. Rotational displacement due to onboard sinusoidal-disturbance moment 

with = 5 and £ = 1 23 

(O n 

12. Rotational displacement due to onboard sinusoidal-disturbance moment 

with — = 1 and £ = 0.25 24 

(O n 

13. Translational single-degree-of-freedom truth model, with base motion 25 


IV 


LIST OF FIGURES (Continued) 


14. Translational displacement due to off-board sinusoidal position disturbance 
with = 0.1 and £ = 0.707 


15. Translational displacement due to off-board sinusoidal position disturbance 

with — = 1 and £ = 0.25 

(O n 


16. Translational displacement due to off-board sinusoidal position disturbance 

with — = 10 and £ = 0.707 
(O n 


17. Rotational displacement due to off-board sinusoidal rotation disturbance 
with = 10 and £ = 0.707 


18. Rotational displacement due to off-board sinusoidal rotation disturbance 
with = 0.1 and £ = 0.707 


19. Rotational displacement due to off-board sinusoidal rotation disturbance 

with = 1 and £ = 0.25 

co n 


25 


26 


26 


27 


28 


28 


v 


LIST OF ACRONYMS 


ARIS 

active rack isolation system 

CAD 

computer-aided design 

IS PR 

international standard payload rack 

ISS 

International Space Station 

TM 

Technical Memorandum 


vi 


NOMENCLATURE 


vector of arbitrary quantity x 

first time derivative of arbitrary quantity x 

unit length of arbitrary quantity x 

arbitrary reference frame x 

arbitrary dynamical system x 

linearized arbitrary quantity x 

arbitrary rigid-body x center of mass 
transpose of matrix x 

ordinary (or total) derivative with respect to time 

partial derivative with respect to scalar u 
trace of arbitrary matrix A 

Scalars 

intersection point on zth actuator arm (fig. 2) 
point locating zth upper stinger (fig. 2) 

point locating zth Lorentz coil (fig. 2) 
system damping matrix: lbf s/ft 


system disturbance input matrix 


NOMENCLATURE (Continued) 


Fi 
F r 
F* 

1 y 

Fi- 


lth 


rAIA' 

l jk 


K 

M 

O 

Q 


Qu 

Qr 



Si 

S u 


magnitude of ith Lorentz coil force vector F C ' : lbf 
point locating ith cross-flexure (fig. 2) 

system holonomic generalized active force, for the rth generalized speed: ftlbf 
system holonomic generalized inertia force, for the rth generalized speed: ft lbf 
system nonholonomic generalized active force, for the rth generalized speed: ft lbf 

system nonholonomic generalized inertia force, for the rth generalized speed: ft lbf 
umbilical attachment point at the flotor end 
stator-fixed reference position of F u 
identity matrix 

central inertia scalar of rigid body A, for body-fixed coordinate directions a_j 
and a k \ ftlbfs 2 

system stiffness matrix: lbf/ft 

system mass matrix: slug, lbf s 2 /ft 

zero matrix 

rotation matrix from the sj coordinate system to the F- ] coordinate system 
ijth element of rotation matrix Q 

contribution to the system holonomic generalized active force for the rth general- 
ized speed, due to rigid body B: ft lbf 

contribution to the system holonomic generalized inertia force for the rth general- 
ized speed, due to rigid body B: ft lbf 

point locating ith lower stinger (fig. 2) 

umbilical attachment point at the stator end 

viii 


NOMENCLATURE (Continued) 


S 


^Fw ^ Fw ^ Fu 

Lowercase 


fj 

fjk 

Si 

k i 

L J 

m A 

Pi 




f jk 


x 0’ 3^0’ z 0 


total dynamical system (stator, flotor, actuators, umbilical) 
geometric lengths: ft 

geometric length for zth actuator assembly: ft 
umbilical damping in the direction: lbf s/ft 
cosine of angle qj for zth actuator assembly 
geometric length for zth actuator assembly: ft 
rotation matrix element 

direction cosines for n ^ in the F\ coordinate system 
umbilical stiffness in the direction: lbf/ft 

geometric length for zth actuator assembly: ft 
mass of arbitrary rigid body A: slug, lbf s 2 /ft 

geometric length for zth actuator assembly: ft 
jth generalized coordinate for zth actuator: rad 
rotation matrix element 

sine of angle z/ / for zth actuator assembly 
jth generalized speed for zth actuator: rad/s 

geometric length for zth actuator assembly: ft 
geometric lengths: ft 


IX 


NOMENCLATURE (Continued) 


x i 

umbilical elongation in the direction: ft 

x Fw yFw Z Fu 

geometric lengths: ft 

x Sw y Sw % Su 

geometric lengths: ft 

K / 

torsional umbilical stiffness about the axis: ft lbf/rad 

< P 

flotor angle of twist relative to stator: rad 

<Pi 

angle-of-twist component in the direction: rad 

Yi 

torsional umbilical damping about the £ axis: ft-lbf- s/rad 

Vectors 


Uppercase 


It 

umbilical bias force in the reference position: lbf 

F C ‘ 

force exerted on the zth Lorentz coil by the flotor: lbf 

F° 

unknown disturbance force acting directly on the flotor: lbf 

F F ' 

force exerted on the flotor by the zth actuator arm: lbf 

F u 

force exerted on the flotor by the umbilical: lbf 

j^AJA* 

angular momentum of arbitrary rigid body A with respect to its mass center A*: 
ft- lbf- s 

m a < 

moment exerted on the zth actuator arm through the upper stinger, due to the zth 
pushrod: ft-lbf 

Mb 

umbilical bias moment in the reference position: ft-lbf 

M c i 

Q. 

moment exerted on the zth Lorentz coil by the flotor, with F ' assumed to act 
at the zth cross flexure: ft-lbf 


x 


NOMENCLATURE (Continued) 


M d unknown disturbance moment acting directly on the flotor, with F D assumed to act 

at the flotor mass center: ft-lbf 

F F 

M 1 moment exerted on the flotor by the ith actuator, with F ' assumed to act at the /th 

cross-flexure: ft-lbf 


Lowercase 


A a B 


d 


i 



C 1 

r AB 


U 

U 1 

AyB 

A a B 


A 


CO 


B 


acceleration of arbitrary point B, with arbitrary reference frame A assumed 
fixed: ft/s 2 

translational acceleration of stator (due to g-jitter): ft/s 2 
system disturbance vector: lbf and ft-lbf elements 
system control current vector: amp 

unit vector in direction of rotation axis for stator-to-flotor rotation (j) 

Vector of generalized coordinates: rad 

position vector from arbitrary point A to arbitrary point B\ ft 

vector of generalized speeds: rad/s 

vector of independent generalized speeds: rad/s 

velocity of arbitrary point B, with arbitrary reference frame A assumed fixed: ft/s 

angular acceleration of arbitrary reference frame B with respect to arbitrary refer- 
ence frame A: rad/s 2 

angular velocity of arbitrary reference frame B with respect to arbitrary reference 
frame A: rad/s 


xi 



TECHNICAL MEMORANDUM 


A “KANE’S DYNAMICS” MODEL 
FOR THE ACTIVE RACK ISOLATION SYSTEM 

PART THREE: ADDITION OF UMBILICALS 
TO THE NONLINEAR MODEL 


1. INTRODUCTION 


As the only active rack-level isolation system for the International Space Station (ISS), the active 
rack isolation system (ARIS) is the central component of the stationwide strategy to meet ISS isola- 
tion requirements. 1 It serves to attenuate the various disturbances that unavoidably accompany manned 
space flight. Umbilicals, as shown in figure 1, are to be used in support of many experiments planned 
for space, providing services such as cooling, power, vacuum, and data transmission. ARIS is nominally 
equipped with thirteen umbilicals. 1 


Vacuum Resource 
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Fire Detection/Maintenance 
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y— Moderate Temperature 
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Figure 1. Detailed diagram of the umbilical assembly. 


Unfortunately, the ARIS umbilicals are not only the primary transmitters of indirect; i.e., off- 
board, disturbances from ISS to the international standard payload rack (ISPR), they are also typically 
nonlinear, hysteretic, and poorly characterized. 2 Since the respective umbilicals are attached at different 
locations to ISS at the bottom plate and to the base of the ISPR at the z-panel, it is desirable to include 
this attachment-point information as part of the ARIS model to improve model fidelity and reduce the 
amount of uncertainty for which the controller must compensate. 


In part one of this series, NASA/TM— 2001-211063, a high-fidelity, linearized, analytical model 
of ARIS was derived using traditional, hand-calculation methods. 3 The model was developed directly, 
using Thomas Kane’s method, without intermediate development of the full nonlinear model. 4 Part one 
briefly outlined an approach for incorporating an ARIS umbilical into the model, using diagonalized 
stiffness and damping matrices. 

Part two, NASA/TM — 2004-213552, presented four computer-based, numerical models of 
ARIS, one of which was purely kinematical, and the remaining three, dynamical. 5 These numerical 
models were used collectively to verify and simplify the linearized analytical model developed in 
part one. 

This Technical Memorandum (TM) completes the nonlinear ARIS model of part two, by adding 
multiple Hookean umbilicals with full stiffness and damping matrices. Comparisons against simple 
single-degree-of-freedom truth models indicate that the completed nonlinear model has input responses 
that follow those expected given the laws of physics. 

This TM (1) briefly summarizes the existing Kane’s model, (2) describes the addition of an 
arbitrary number of massless umbilicals with specified arbitrary attachment points and arbitrary (paral- 
lel) stiffness and dampings, to that model, (3) provides a brief description of Autolev™ in which the 
underlying untethered model was developed, and (4) presents the process by which verification of the 
enhanced model was accomplished along with some verification results. 

To develop the enhanced model, it was assumed that the umbilicals were to be included individu- 
ally in the model, rather than as a combined, or effective, umbilical. For the sake of simplicity and since 
the range of motion is small, each umbilical is assumed to be Hookean and massless since an algebraic 
state-space model, one that is useful for controller design using state-space-based methods, is desired. 
Also, each is modeled as a three-dimensional (six-by-six, constant) stiffness matrix in parallel with a 
three-dimensional (six-by-six, constant) damping matrix. 

In the verification stage, the model of the tethered ARIS system is compared with two baseline, 
or truth, models. One truth model consists of a translational, one-dimensional, second-order spring- 
mass-damper system, and the other consists of a corresponding rotational system. 6 The verification 
checks indicate that the proposed Kane’s model, when properly calibrated against experimental data, 
will provide a high-fidelity representation of ARIS for simulation and controller-design purposes. 
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2. COORDINATE SYSTEMS 


The coordinate systems, generalized coordinates, generalized speeds, and angular velocities used 
for the Kane’s model were presented previously in part one. 3 ’ 7 ’ 8 They are reviewed in this section and in 
sections 3 and 4 for the reader’s convenience. 

With the IS PR (flotor) in the reference position; i.e., centered in its ISS- fi xcd rattle space, fix 
eight right-handed, orthogonal coordinate systems in the ISPR, one at each of the cross-flexure cen- 
ters (fig. 2). Let the zth coordinate system have origin F { , located at the center of the zth cross flexure, 

V 

(z=l,...,8), with axis directions determined by an orthonormal set of unit vectors, J_j (/’=1 ,2,3). The 

overhat indicates unit length, the index, z, corresponds to the zth actuator assembly, and the index, j, 

2-i 

distinguishes the three vectors. Orient the unit vectors such that / 9 is along the zth arm, toward the zth 

p 

voice coil; / j is directed parallel to the other segment of the zth arm and toward the upper stinger, which 
is located at A 2 ; and J 3 is in the direction /j x / 2 along the intersection of the two crosspieces of the 
zth cross flexure. 


Coil 



Figure 2. Kinematic diagram, including the zth actuator assembly and the umbilical. 
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Fix a similar right-handed coordinate system, a l -(j= 1,2,3), in the arm of each actuator. Locate 
each system, a!-, such that it is coincident with the corresponding flotor-fixed coordinate system, f_j , 
when the flotor is in the reference position. 


Locate eight more arm-fixed coordinate systems, described by a’j , at the respective upper sting- 
ers, points A 2 . At the respective lower stingers, points S-, place eight pushrod-fixed coordinate systems 

p\ , and eight stator-fixed (/SS- fixed) coordinate systems, s 1 - . Orient these 24 coordinate systems such 
— J J 

that when the stingers are relaxed; i.e., with the flotor in the reference position, the coordinate directions, 
a 1 - , p \ , and s'-, are coaligned for the /th actuator, with p 1 (along with a' 2 and s' 2 , in the reference 
position) directed from 5) towards A 2 . Directions p* and p 1 ^ (in the reference position) are determined 
by fixed rotation matrices relating the a - to the cx - coordinate systems. 

Finally, define a primary, central, flotor-fixed, reference coordinate system with coordinate direc- 
tions, f_j . All other flotor-fixed coordinate systems, / ■ , are assumed capable of being referenced; e.g., 
by known direction cosine angles, to this system. 


Define a stator-fixed coordinate system, , that is associated with the jth umbilical. Define the 

~~ l ~ j 

rotation matrix from the stator- fixed coordinate system, , to the new coordinate system, , by 



( 1 ) 


where 



fit 

>\2 

r 13 

J 

r 2l 

r 22 

r 23 


r 32 

r 33 


( 2 ) 


Note that |^9 Xj J 


is a constant matrix, and it could differ for different umbilical circumstances. 
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3. GENERALIZED COORDINATES AND GENERALIZED SPEEDS 


Let the a - coordinate system rotate, relative to the / ■ coordinate system, through the positive 
J V ' 

angle, q{, about the / 3 axis against the cross-flexure stiffness, k \ . Similarly, let the orientation of the 
a'j coordinate system, relative to the p 1 . coordinate system, be described by consecutive positive rota- 
tions— q\ about the p l axis and q\ about the moved 3 axis. Also, let the orientation of the p l . coordi- 
nate system, relative to the s'j coordinate system, be described by consecutive positive rotations — q\ 
about the s\ axis, q l 5 about the moved 2 axis, and q l 6 about the moved 1 axis. The six generalized 
speeds are defined as the time rates of change of the respective generalized coordinates: 

u ) =<7y(/ = l,-",8;j' = l,--,6) . (3) 

Finally, let c l j and s l j represent the cosines and sines of the respective angles, q'j . 
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4. ANGULAR VELOCITIES OF REFERENCE FRAMES AND RIGID BODIES 


Designate the reference frames corresponding to the stator, the zth pushrod, the zth arm, and the 
flotor, as the symbols S, P n A,-, and F , respectively. Let 5, and F • represent, respectively, the coordinate 

. -iT 


systems in S and F, respectively defined by 


. nT 


*1 -2 


and 


L U L 


. Two intermediate 


reference frames were introduced previously to permit describing the angular velocity of each pushrod 
relative to the stator. Designate those intermediate frames corresponding to the zth actuator assembly as 
Rj and Q r Another intermediate reference frame was previously introduced between frames P t and A,-; 
designate that frame as 7) . Finally, designate the stator-fixed reference frame corresponding to the /th 


umbilical as Ej . 


Let each intermediate reference frame have a frame-fixed, dextral set of unit vectors. Indicate 
the unit vectors for each of these frame-fixed coordinate systems by using the corresponding lowercase 
letter; i.e., P- corresponding to R n etc. The following gives the expressions for the angular velocities 
of the various reference frames and rigid bodies of S : 


F 'C0 Ai =u[l^ f 

(4) 

P <0) T ' =u l 2 P x , 

(5) 

T A i *i 

'r‘ = «3l 3 . 

(6) 

s i 0 ) R i = u\ S 3 , 

(7) 

R ‘ ofi‘ =U l 5 ?2 , 

(8) 

Q '0) P ' =U l 6 q\ . 

(9) 


Using the addition theorem for angular velocities, the angular velocities of the rigid bodies 
of S are 


S ■ A ■ / /v / / i M i /V/' / a / 

'CO ' = Uo p + j L ^4 ^3 +zz^r^ -I - u^C[^ 


( 10 ) 
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S '• --P: i '•i , i '•i . i '•i 

1 (0 1 — M 4 s_2 + M 5 + UfoC] , 


( 11 ) 


and 


F. / a / 7 A/ / a / 

1 OJ 1 = 1/4 s ^ + M5 + Mg <7^ • 


( 12 ) 
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5. GENERALIZED ACTIVE FORCE CONTRIBUTIONS DUE TO THE UMBILICALS 


5.1 Generalized Active Force Equations 

Equations of motion developed using Kane’s method consist of contributions called generalized 
active forces (GAFs), related to system forces and moments, and contributions called generalized inertia 
forces (GIFs), related to the time derivatives of linear and angular momenta . 4 For ARIS, the develop- 
ment of the GAF and GIF contributions due to the actuator arms and pushrods, the flotor, and a single 
umbilical with diagonal stiffness and damping matrices is detailed in part one . 3 Because the umbilicals 
are assumed to be massless, they make no contributions to the GIFs, and thus, they contribute only to 
the GAFs. 

In previous treatments (part one of this series and a paper by Rupert and Hampton), the umbilical 
stiffness and damping matrices were assumed to be diagonal with no coupling between translation and 
rotation, or between axes . 3 ’ 8 However, coupling will always exist between translation and rotation and 
typically between axes as well . 3,9 In other words, for an arbitrary set of orthonormal axes, one expects 
that the six-by-six umbilical stiffness and damping matrices will be full. Consequently, full matrices are 
assumed in the following treatment. 

Designate the rth partial velocity of the /th umbilical’s flotor- attachment point (with reference 
frame S assumed to be fixed, for purposes of differentiation) as 1 v u i and the rth partial angular 

c zp 

velocity of the flotor relative to the stator as co r . Designate the force and moment that the jth umbili- 
cal exerts on the flotor as F U J and M U J, respectively, where the force is assumed to act on the flotor at 
the jth umbilical attachment point F V j. Fet all direct disturbance forces and moments that are exerted 
on the flotor— such as those due to air currents, to direct contact by astronauts, or to moving parts on 
flotor-mounted experiments— be designated collectively by the respective symbols FP and MP ■ Finally, 
let - Fp ' and -M Ci represent, respectively, the force and moment exerted by the /th Forentz coil, located 
at A 3 , on the flotor, where 


A' 


• C 'q\ 


(13) 


assumed to act at point F i and 


M Ci = r_ FiA 3 x F Ci = -F C ' (l[ + l\ )a ' 3 . (14) 

Then, assuming a nominal complement of 13 umbilicals, the GAF contribution due to the flotor, for 
the rth generalized speed, can be written as 


8 


( 15 ) 


13 


Qr := Sr' -r Uj ' ^ Jj + V ' K U 
j= A 

8 


+ 5 «f 


M D + '^ J k[q\[ 


l 

—3 


i = 1 


+ s i v f 170 


-? ‘ - C ' + 5l ®f ■ M c ' ) • 

The /th umbilical affects only the first two terms, with the contribution 

7 Q r F := v_ F r u i ■ F U J + s a) F ■ M_ U J . 
These two terms will be treated in sections 5.2 and 5.3. 


(16) 


5.2 Partial Velocities Sl v Fu J 


Represent by r the position vector from arbitrary point A to arbitrary point B, and define 
the following position vectors using the indicated scalars: 


, F i A 2 _ li ;.i , 7 i ~ ' 

— ci 2 t“ ? 

(17) 

SjAn ,i ~ i 

r ‘ =hP 2 > 

(18) 

i ~i 

r =P2P 2 . 

(19) 

UA; i * i i ~ i 

1 1 = ^ 

(20) 


and 


FF 


l 


/ 1 / 1 +/ 2 / 2 +/ 3/3 • 


( 21 ) 


Represent the position vector from the flotor center of mass to the /th umbilical’s flotor- attachment 
point as 

zp* zp i * 1 / 1 i * 1 

r b b Ui -yi f xW f x7i " 


r u j 


xi /, +r 7 +z 7 

F u J -\ f u j -2 FjjL 3 


( 22 ) 


for appropriately defined measure numbers. Then, referring to figure 2, the position vector from 
the first-lower-stinger attachment point to the flotor center of mass is seen to be 


= i\p\ - /tea - + fi l f\ + f\? 2 + fif\ • 


(23) 
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Since 


r S i F Uj = r' S|/ + r F F Uj , 


the rth partial velocity of the flotor-attachment point F U ; is 6 


c Fjj . 

S i v J 


Su.F 


>v + 


a ,s > /i x A 


5m. 


(24) 


(25) 


c c 1 " F F S F 

The position vectors r 1 and r u j, and the angular velocity vector 1 m 1 , can be expressed 

in terms of the coordinates, qj (j= 1,...,6), and the associated generalized speeds of actuator number 1. 
This means that the velocity of each attachment point, F u ., can be expressed in terms of these same 
coordinates and generalized speeds so that the partial velocities defined in equation (25) are nonzero 
for only those six generalized speeds; i.e., for r=l,...,6. 


£ J 

Expanding equation (25) in terms of the /th stator-fixed coordinate system, <§ . , (associated with 
the jth umbilical) and collecting terms, the rth partial velocity for can be expressed as 


Si 



3 



i = 1 


(26) 


where the scalars, v ri j , are simply the measure numbers in analytical (algebraic) form. In matrix form, 
the six nonzero partial velocities for the /th umbilical’s flotor-attachment point are 


Si 

Si 

Si 

< 

Si 

Si 

Si 



i=t 

3 

X u, y 

V 2 i'l 

i = 1 

i xn 

i = 1 

i= t 

i= 1 
3 

X Uj y 


J= 1 


Uj 

V 11 

Uj 

V 12 

U j 
v 13 

Uj 

V 21 

Uj 

v 22 

t 7 / 

v 23 

Uj 

V 31 

Uj 

v 32 

U J 

v 33 

Uj 

V 41 

Uj 

v 42 

t 7 / 

v 43 

Uj 

V 51 

t// 

v 52 

U J 

v 53 

Uj 

V 61 

Uj 

V 62 

t 7 / 

v 63 


^ / 
S 

l’ 

—2 


-3 


(27) 
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5.3 Partial Angular Velocities S ' 0 J r Fl 

From equation (12), the angular velocity of the flotor relative to the stator can be expressed as 


S, F, 1 Id 1-4 1 -4 l l rt 

(0 — u 3 p j + ^ + u 4 £3 + W 5 ~ U\ f 3 


( 28 ) 


Following the procedure and using a notation analogous to that of the previous subsection, the nonzero 
partial angular velocities can be expressed in the matrix form, 



Uj 

®n 

U j 

°>12 

1 

3 

Uj 

®21 

Uj 

( 0 2 2 

Uj 

( 0 23 

Uj 

®31 

Uj 

®32 

Uj 

(0 33 

Uj 

co 4 { 

Uj 

®42 

Uj 

®43 

Uj 

®51 

Uj 

®52 

Uj 

®53 

Uj 

®61 

Uj 

®62 

3 

22= 

1 


~ / 
i 

* j 
-2 

* j 
-3 


(29) 


5.4 General Form for Umbilical Forces and Moments 


With full stiffness and damping matrices, the /th-umbilical force, F U J, and moment, M U i, will 
each be expressed in terms of umbilical-elongation-, elongation-rate-, angle-of-twist-, and twist-rate- 
component measure numbers. In vector form, the umbilical elongation is r Uh j 11 j, where F uh desig- 
nates the stator-fixed location of the umbilical attachment point, F f j, in the reference position. Let (pn^ 
represent, in vector form, the rotation of the flotor, relative to the stator, from the reference position. 

The term n ^ is in the positive direction of the rotation axis, and 0 is the positive angle of twist about 
that axis. Then, for the jth umbilical, the elongations; i.e., the elongation measure numbers, are 

x{ =r Fuh j Fu j -I 7 , (30) 

1 ~ —l 

the elongation rates are xj , the angles of twist are 


11 


( 31 ) 



and the twist rates are 0/. 

Let the bias force and moment exerted on the flotor by the /th umbilical be F_ b 1 and M_ b 1 , 
respectively. Next, for z'=l,2,3, define 




M bi = M^ j ’I- . 

& 




and 


M, , fUi M 


M: 1 = M u j ■ £ . 

' — — i 

Then, in matrix form, the combined force and moment equations are as follows: 


M± j 

M^ j 


k Uj 

K \ l 


.61 


k U J 

k \6 


k U J 

k 66 


K 


U J 

x \ 

u i 
x 2 1 

U, 


<t>l j 

$ 


1 + 


Uj 

11 


Uj 

-61 


Uj 

16 


C U J 

C 66 


c 


' .Ui' 


f Ui ) 

X 1 


hi 

■ Ui 


Ui 

x 2 


hi 

.U: 


Ui 

x 3 


hi 

> 

■Ui 


Ui 

01 7 


M b { 

■Ui 


u i 

<P 2 ] 


M b2 

■Ui 


U i 

w\ 


Wbi \ 


(32) 

(33) 

(34) 


(35) 


(36) 


where the stiffness and damping elements, obtained by measurements or other means, and the stiffness 
and damping matrices are indicated by the designations k, c, K, and C, respectively. The stiffness and 
damping element superscripts and subscripts have their obvious meanings. To use the left-hand elements 
of equation (36) in equation (16), these elongations and angles of twist, along with their rates, must now 
be reexpressed in terms of the system coordinates and generalized speeds. 


5.5 Umbilical Elongations and Elongation Rates 


Recall equation (30), for z'=l,2,3; but 

r F uhj F Uj =L S Mi +r _ F ' F Uj -r S i s uj -r S Uj F uhj j 


(37) 
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where the right-hand-side terms can be expressed by 


S,F t A -1 ,1 -1 ,1 -1 

r 1 1 — I 3 P 2 ~ h(±i ~~ h ’ 

(38) 


(39) 

r ~ s,Sv i = 4 „ + >i i 2 + 4, i'i • 

(40) 


and 


r Su J Fuh J =4l\ + y J 0^2 + Zol\ 


(41) 


for appropriately defined coefficients. 

~ j 

Using equation (1), equation (37) can now be written in terms of the £ . coordinate system in lin- 
earized form 


r F Uhj F Uj - r J 


= x{^ + x^_ J 2 + x^_ 


(42) 


where 


j 

X 1 


c{ + y J Fu (<7l ~ <73 ~ 74 ) + Zp v <?5 A (73 + <?4 ) 

J 

x 2 

II 

c i - (^2 + <?6 ) _ X F V ( 9 ! _ <?3 - <74 ) ~ 4 (<?3 + c l\ ) 

J 

x 3 


C{ + [l\ ~ + |y4 - A )(<72 + <76 ) + A<l6 


(43) 


for 


and 


i— < cs 

K 

1 

K 

II 

" ——a 

0 

(44) 

C 2=y J F u - >’4 + ; 3 - A - >’() . 

(45) 

c 3=4 u -4 u - z o ■ 

(46) 


Differentiating equation (43)/ 
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( 47 ) 


V 

X 1 


y J Fu - u\ - u\ ) + z J Fu u\ -l\u\ + 1\ (u\ + u \ ) 

x j 

x 2 

II 

~ Z F U ( u 2 + u 6 ) - X F V ( M 1 - u 3 - u 4 ) - 4 [ u 3 + u 4 ) 

x 3 


[ l 2 ~ X F U )“5 + [y J F u -^l)( M 2 + m 6 ) + ^3 U 6 


Equations (43) and (47) give the umbilical elongations and elongation rates in terms of the six indepen- 
dent generalized speeds. It is now required to do the same for the umbilical angles and angles of twist. 

5.6 Umbilical Angles of Twist and Twist Rates 


Express h ^ as 


Define rotation matrix Q by 


✓s a 1 a 1 a 1 

^ = Slh+g2l2+S3l3 • 


/' 


f s 1 j 

S-I 

/' 

—2 

II 

jo 

^ 1 

*2 

f 

—3 


S 1 

-3 


(48) 


(49) 


and then the linearized three-by-three rotation matrix, t Q, has elements, jQ^, defined as 


f 

J -l 

r 

* 1 


U 

" = 

* 1 


[/ 3 J 

- 


4i 


q\ 

45 


1 , 1,1 

"4 1 +43+44 


44 


“42 


46 


1 

-45 
42 +46 


“ 

f -1 ] 


h 


* 1 

< 

-2 


^ 1 


[-3J 


(50) 


Let the postsuperscript T indicate matrix transposition and let tr /Q represent the trace of jQ. For small 
angles it can be shown that 10 


0 

S3 

S2 

S3 

0 

8l 

§2 

Si 

0 


iQ~iQ r 

(i+tr/e) 1/2 


( 51 ) 
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Substitution for / Q , as defined in equation (50), into equation (51) and simplification yield 


and 



(52) 

(53) 


1 



(<7l “ c l\ ~ c l \ ) • 


(54) 


Substituting from equations (52)-(54) into equation (48) and transforming into the sj coordinate sys- 
tems by use of /Q yield the following linearized expression for the flotor rotation: 


i / \ ( 1 ] \ /v ] 1 -*.] / 1 1 1 \ / C C\ 

^^={^ 2 + q 6 )s_ x + q 5 s_ 2 -\c lx - . (55) 

Use of equations (31) and (48) leads at last to the linearized forms for angular position and rotation rate: 


<p\ 


q\ + < ?6 

02 

-M 

<75 

03 


- (<7i - q\ - q \ ) 


01 


u\ +u\ 

02 


u\ 

03 


— — U 3 — M 4 j 


(56) 


(57) 


The jth umbilical portion of the GAF contribution due to the flotor for the rth generalized speed 
can now be expressed in terms of the coordinates and independent generalized speeds: 
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where 


J Q F = 5 ‘ v, Uj ■ F U J + s co F r ■ M U J 


X 

r = 1 
6 

I 

r = 1 
6 

I 

r=l 

6 

I ' 

r=l 

6 

X < 

r=l 

6 

X < 

,r=l 


U J 

; r\ 


U J 

’A 


Uj 

’A 


u i 

W 


u. 


u. 


& 

F U A 

M^ j 


(58) 


& 

m\' 

m u 2 > 


k u j 


k u J 

.61 


*16 


*66 


*1 

Uj 

x 2 

U, 


4>i J 

$ 

< t> u A 


1 + 


U] 

ii 


u j 

~61 


16 


'66 




' JU. ' 

X 1 


F bl 

.[/, 


1 Ui 

x 2 


F l:2 

■ Ui 


Ui 

x 3 


F b3 

> 

■Ui 


Ui 



M b { 

■Ui 


Ui 



m,A 

■Ui 


Ui 

w\ 


[Mbi \ 


(59) 


and the elongations, elongation rates, angles of twist, and twist rates are as defined in equations (43), 
(47), (56), and (57). The model is now ready for calibration and computer implementation. 
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6. IMPLEMENTATION IN AUTOLEV 


A nonlinear, rigid-body model for the untethered ARIS; i.e., without umbilicals, with full 
actuator dynamics (masses and inertias), was previously developed using OnLine Dynamics’ Autolev, 
a DOS-based interpreter specifically designed to solve dynamics problems using Kane’s method. 11 
The umbilical contributions developed above (eqs. (58) and (59)) were incorporated into that model 
by appropriately adding those contributions to the untethered-model GAFs. In Autolev, this required 
entering the locations of the umbilical attachment points, the orientations of the stator-fixed coordinate 
systems associated with the umbilicals, and the umbilical forces and moments determined above. 

The geometric data required to include the umbilicals were determined from CAD models, such 
as in figure 3. Since umbilical stiffness and damping values were not available, representative estimates 
were entered and varied as needed for the process of model verification (described in section 7). 



Figure 3. CAD-based technique for determining umbilical attachment locations. 


After adding necessary terms to the Autolev model of ARIS, the code was compiled in the 
Autolev workspace. Autolev then auto-generated a C-program file containing the equations of motion 
and a numerical integration routine for their solution. Results of executing the C-program were plotted 
using Microsoft® Excel. 
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7. MODEL VERIFICATION 


Numerous simulations were run to verify the model for special cases. The basic procedure was 
to reduce the actuator stiffness to very small levels and apply selected test forces (moments) through 
(about) the flotor center of mass. Corresponding simulations were compared between the full nonlinear 
model and appropriate (translational or torsional) single-degree-of-freedom truth models. Representative 
comparison checks are described in sections 7. 1-7.8. 

For each of the following comparisons, a single umbilical is assumed to be attached between the 
stator and the flotor center of mass. Umbilical biases are assumed to be zero. The umbilical coordinate 
system, <§ ., is assumed to be aligned with the flotor-fixed coordinate system, ./ . , when the flotor is in the 
reference position. Stiffness and damping terms that couple translation and rotation are all assumed to be 
zero. For forces (moments) applied through (about) the flotor mass center and for small-enough actuator 
stiffness, the ARIS system dynamics should approximate those of a simple second-order spring-mass- 
damper system (fig. 4). 



k Translational Stiffness 
c Translational (Viscous) Damping 
m Mass 
F Applied Force 
x Translational Displacement 

kx ex 


m 

T 

F 


Figure 4. Translational single-degree-of-freedom truth model. 


The flotor mass is assumed to be 55.6 slugs (1,790 lbm, 811.94 kg), and the central principal 
moments of inertia, 166.8 slug-ft 2 (5,371 lbm-ft 2 , 226.34 kg-m 2 ). The umbilical stiffness matrix is 
assumed to be diagonal with each diagonal element set at 30 (lbf/ft for translational elements and 
lbf-ft/rad for rotational elements). The damping matrix is also assumed diagonal with all diagonal 
terms numerically equal but varied as needed to control the damping ratios, depending on the partic- 
ular simulation. 
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7.1 Onboard Impulsive-Disturbance Force, No Damping 

For the first test of the ARIS model, an impulsive disturbance force was applied directly to 
the flotor mass center along the £ direction, with zero damping. The disturbance was also applied 
to a simple spring-mass truth model with the same mass and stiffness. The impulse was approximated 
by a rapidly decaying exponential function: 


F = Ae~ dt i 2 , (60) 

where d= 1 00 s -1 and A=100 lbf (444.82 N). The two systems had identical sinusoidal responses (fig. 5) 
at the expected natural frequency of 0.7346 rad/s (0.1169 Flz). 



Figure 5. Translational displacement due to onboard impulsive-disturbance force. 


7.2 Onboard Sinusoidal-Disturbance Force, With Damping 

In a second test, a sinusoidal disturbance force with an amplitude of 3 lbf (13.34 N) and a fre- 
quency of ©was applied to the flotor mass center, again in the ^ direction, using various frequency 
ratios, co/co n , and damping factors, £. As before, the disturbance was also applied to a corresponding 
truth model. The two system responses were identical (figs. 6-8). 
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« ARIS and Truth Model Results 
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Figure 6. Translational displacement due to onboard sinusoidal-disturbance force 
with = 1 and £ = 1 . 



Time (s) 


Figure 7. Translational displacement due to onboard sinusoidal-disturbance force 

with = 5 and £ = 1 . 

(O n 
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Figure 8. Translational displacement due to onboard sinusoidal-disturbance force 
0) 

with — = 1 and £ = 0.25. 

co n 

7.3 Onboard Impulsive-Disturbance Moment, No Damping 

To check the ARIS model rotational response, an impulsive disturbance moment was applied 
directly to the flotor at its mass center about the direction with zero damping. The disturbance was 
also applied to a simple rotational spring-mass truth model, with the same inertia and rotational stiffness. 
Again, the impulse was approximated by a rapidly decaying exponential function: 

M=Ae~ dt i 2 , (61) 

where d= 1 00 s -1 and A=100 lbf-ft/s (135.58 N-m/s). The two systems had identical sinusoidal responses 
(fig. 9) at the expected natural frequency of 0.4241 rad/s (0.0675 Hz). 
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Figure 9. Rotational displacement due to onboard impulsive-disturbance moment. 


7.4 Onboard Disturbance Moment, With Damping 

To test the damped rotational response, a sinusoidal disturbance moment with amplitude 
0.09 lbf-ft/s (0.122 N-m/s) and frequency, ft), was applied at the flotor mass center, about the direc- 
tion, using various frequency ratios, co/co n , and damping factors, £. The disturbance was also applied 
to a corresponding truth model. Again, the two system responses were identical (figs. 10-12). 
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Figure 10. Rotational displacement due to onboard sinusoidal-disturbance moment 
with = 1 and £ = 1 . 



Time (s) 


Figure 11. Rotational displacement due to onboard sinusoidal-disturbance moment 
with = 5 and £ = 1 . 
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Figure 12. Rotational displacement due to onboard sinusoidal-disturbance moment 

with = 1 and £ = 0.25 . 

(O n 

7.5 Off-Board Translational Disturbance 

In normal operation, the flotor will also be subject to off-board disturbances, which are transmit- 
ted via the ARIS umbilicals and actuators. With negligible actuator stiffness, the system response should 
match that of a corresponding truth model (fig. 13). A translational sinusoidal disturbance of amplitude 
0.25 in (0.635 cm) was applied to the stator (base) in the ^ direction with the ARIS umbilical oriented 
so that its stator attachment point was collinear with the flotor mass center, along the c, 1 direction. (The 

flotor attachment point was assumed to be at the mass center, as before.) Various frequency ratios, co/co n , 
and damping factors, £, were used with the flotor- to- stator relative displacements compared to those of 
the truth model. The responses of the two systems were identical (figs. 14-16). 
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Figure 13. Translational single-degree-of-freedom truth model, with base motion. 
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Figure 14. Translational displacement due to off-board sinusoidal position disturbance 

with = 0.1 and £ = 0.707 . 
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15. Translational displacement due to off-board sinusoidal position disturbance 


with = 1 and £ = 0.25. 



Figure 16. Translational displacement due to off-board sinusoidal position disturbance 


7.6 Off-Board Rotational Disturbance 


A rotational sinusoidal disturbance, of amplitude 0.001 rad (0.0573°) was applied to the stator 
(base) about the £ direction with the ARIS umbilical attached and oriented as in the previous section. 
Various frequency ratios, co/co n , and damping factors, £, were used with the flotor-to-stator relative rota- 
tions compared to those of the truth model. The responses of the two systems were identical (figs. 17-19). 



Time (s) 

Figure 17. Rotational displacement due to off-board sinusoidal rotation disturbance 
with = 10 and £ = 0.707. 
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Figure 18. Rotational displacement due to off-board sinusoidal rotation disturbance 

with = 0.1 and £ = 0.707. 
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Figure 19. Rotational displacement due to off-board sinusoidal rotation disturbance 
with (0 - = 1 and £ = 0.25. 

co n 
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7.7 Comparisons Using Nondiagonal Stiffness and Damping Matrices 


Suppose now that the stiffness and damping matrices, [ A'] and [C] , respectively, used for the 


above simulations are changed by selecting a new coordinate system, rotated through 
the previous coordinate system, | : 


91 


from 


i; 

I - i — 

1, 

L 

II 

«2 

1 

-3 


l 3 . 




(62) 


Then the upper left and lower right block diagonal three-by-three submatrices, [^i i ] and \Kj2 ] , of the 
(initially diagonal) stiffness matrix, [/f] , become, respectively, 


and 





(63) 





(64) 


Similar relationships hold for the damping matrix, [C] , and its corresponding submatrices. If the 
same disturbance force or moment is now applied as in any of the preceding simulations, then the result- 
ing system motion should remain unmodified, which was shown to be the case . 12 This indicates that the 
ARIS model correctly handles cross-coupling among axes; i.e., full stiffness and damping submatrices, 

[KllM x 22].[ c ll].and[C 22 ]. 
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8. CONCLUSION 


This TM has shown how multiple ARIS umbilicals, modeled as massless three-dimensional 
Hookean springs in parallel with three-dimensional viscous dampers, are to be included in the nonlinear 
Kane’s model of ARIS. The umbilical force, moment, partial-velocity, and partial-angular- velocity terms 
were determined for the flotor contributions to the GAFs in Kane’s dynamical equations. 

Next, these modifications were incorporated into the ARIS Autolev simulation code. Simplify- 
ing assumptions were made to permit the comparison of the Autolev model’s input responses to those 
of simple second-order truth models for selected test inputs. Numerous simulations, involving onboard 
and off-board force and moment disturbances, showed that the ARIS-model responses match those of 
the truth models. This indicates that the proposed model, when properly calibrated against experimental 
data, will provide a high-fidelity representation of ARIS for simulation and controller-design purposes. 
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